FANSE-96/09 
LPTHE Orsay-96/24 
lhep^ph/96042731 

Relativistic Quantum Scattering of High Energy Fermions in 
the Presence of Phase Transition 

Jose Rodrfguez-Quintero and Manuel Lozano 

Departamento de Fisica Atomica, Molecular y Nuclear, Universidad de Sevilla, SpainQ 

Olivier Pene 

LPTHE, F 91405 Orsay, France^ 

Abstract 

We study the high energy behaviour of fermions hitting a general wall caused 
by a first-order phase transition. The wall profile is introduced through a general 
analytic function. The reflection coefficient is computed in the high energy limit 
and expressed in terms of the poles of the wall profile function. It is shown that the 
leading singularity gives the high energy behaviour. 
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Much work have been devoted to the problem of transmission and reflection of rela- 
tivistic fermions through a wall separating two phases of different symmetry properties. 
The main effort concentrates in developing the idea |jl|] that the baryon asymmetry of the 
Universe might have been produced if the cosmological electroweak {SU{2) x U{1)) phase 
transition has been of first order. In these works the phase transition is described in terms 
of bubbles of "true" vacuum with an inner expectation value of the Higgs field f 7^ 0, 
i.e. a spontaneously broken symmetry phase, appearing and expanding in the preexisting 
"false" vacuum with v = 0, i.e. a symmetric phase. 

In this scenario the quarks/ ant iquarks hitting the wall from the unbroken phase are 
reflected or transmitted. The point to elucidate is whether there exists a CP-asymmetry 
that produces a different reflection and transmission probability for quarks and antiquarks 
in order to explain, via the standard model baryon number anomaly 0, 

d,r^ = d,r, = Nf{^ww - ^YY). (1) 

the correct baryon asymmetry of the Universe. In the physical conditions of the early 
universe the fermions moving through the bubble wall will interact not only with the 
wall but also with the particles in the surrounding plasma, thus we have a transport 
problem. This transport plays an essential in the so-called "charge transport" mechanism 
0], in which the action of the baryon anomaly happens at a distance from the bubble 
wall. This diffusion problem is very complicated and involves solving the Fokker-Planck 
equation taking into account CP violation and baryon anomaly. A useful simplifying 
assumption is to decompose the process into two steps, one describing the production 
of the CP asymmetry when the quarks/antiquarks are reflected on the wall, the second 
describing the transport and the eventual transformation of the CP asymmetry into a 
baryon asymmetry via the baryon number anomaly (|l|). 

One expects that the diffusion corrections are relatively minor ones to the first step, 
i.e. to the scattering from the wall, although effects of the surrounding plasma are incor- 
porated by mean of introducing the Higgs field effective potential that takes into account 
the temperature of a thermal bath. The structure of the wall depends on this effective 
potential and its knowledge is obviously necessary in order to compute the reflection and 
transmission coefficients. However, the profile obtained by solving the equation of motion 
is rather complex and depends on many coupling constants [^. At this point, two lines 
for simplifying the problem have been followed. The first is to describe the profile of the 
bubble wall by an analytical function that simulates the dynamics of the phase transition 
and then to treat the scattering in an approximate way [1-3]. The second approximates 
the wall profile by a step function i.e. a sudden jump from one phase to the other, the 
thin wall approximation 0, [0. this extreme simplification allows to compute in an exact 
way the Feynman fermion propagator in the presence of a wall 0. 

Our aim in the present paper is to study the general problem of fermions hitting a 
wall in order to connect the analytic properties of the profile function to the behaviour 
of the fermions. In general, this question is very difficult, but interesting conclusions 
about this relation between the proffie and the behaviour of fermions in the high energy 
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limit can be obtained. We develop a general method to calculate the reflection and 
transmission coefficients of high energy fermions hitting a wall, establishing a relationship 
between the poles of the profile function and these coefficients. The quantum corrections 
to the expected classical behaviour are obtained. Apart from the clear theoretical and 
purely formal interest of this result, cosmological imphcations and possible applications 
of the formalism to systems of relativistic fermions undergoing a phase transition should 
be considered. Some implications in CP- violating process are studied in Q. Two other 
examples of application fields can be condensed matter under extreme external conditions 
or certain stages of the quark-gluon plasma formation process. 

As usual in this kind of work, we formulate the problem in the rest frame of a wall, 
parallel to the x — y plane and normal to the z-axis, characterized by a general profile. 
In order to calculate the reflection coefficient we need only the plane wave solutions of 
the Dirac equation for particles moving along the z-axis. In any case, for other incoming 
directions we can perform an appropriate Lorentz boost parallel to the x — y plane and 
reduce the problem to the latter. The phase transition is incorporated into the Dirac 
equation by including a position dependent mass term which varies inside a certain region 
designated as the domain wall and takes two different constant values in the two outer 
sides of the wall. Following Nelson et al[0, we work in the chiral basis, conveniently 

reordered to obtain 7^ = f ^ \ and factor the Dirac operator into 2x2 blocks. 



-ag. 

Thus the Dirac equation can be expressed as 

(id, + idt -m{z) ^ \ 

m* (z) idz — idt 

idz + idt —m*{z 

V m{z) idz — idt J 

With the following Ansatz for solutions with positive energy E 



^ = 0. (2) 



^=('^']e-'' with ^^-(|) (3) 

where ipi and '?/'4 are eigenspinors of the chirality operator, 75, for the eigenvalue -|-1 and 
iIj2 and ips for —1, we obtain 

itdz + Qiz))^j = 

{tdz + Q{z))^n = (4) 

being 

<2W = f,„L "'tM -d <3W = r„l -'"W'l; (5) 



where the mass function can be considered as complex in order to incorporate CP- violating 
process in the formalism. In this work we are not interested in CP-violation but in a 
general approach of the problem. Thus the mass function will be assumed as real and 
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consequently Q{z) and Q{z) will be identical matrices. The solution for the first equation 
of can be written as follows 

Mz) = Pe ^^(^^^ j (6) 

and analogously for the second one 

v,„(.)^p.-^;-^«(:{;:j:j), (7) 

where P indicates a path ordered product and r is the position variable along the z-axis. 
Nevertheless, it is obvious that for a real mass function the quantity 

SJ(2,2„) = Pe'J'-o**" (8) 

in © i= the =ame than = Pe^;-^«W in ©. In consequence, the task i. to 

evaluate equation (H). 

Using the usual Pauli's matrices, Q{t) can be expressed as 

Q{t) =(T3E -ia2m{T). (9) 

We consider now 



m(r)=mo/(r), (10) 

where /(r) is a certain function which describes the structure of the domain wall, the 
profile wall function. The asymptotic conditions /(+oo) = 1 and /(— oo) = are required, 
/(r) — 9{t) decreasing exponentially when r ±oo. it is also assumed that /(r) = 0(1). 
The profile function will be considered as an analytic function in the real axis, therefore 
a profile which is constant outside a certain finite region, the domain wall, cannot be 
described. This last kind of wall profiles are studied in detail in If the profile function 
is analytic, the wall, defined as the region where the mass varies with the position, extends 
formally from r = — oo to r = +oo. Nevertheless the domain wall will be characterized 
by an effective thickness, a. The particular criterion considered in order to define this 
parameter is not important; for instance, it can be established by taking into account that 
|/(r) — 1| < 0.1 for T > a and |/(t)| < 0.1 for r < —a. However, in general we can write 

/M = i^(^) • (11) 

A characteristic energy, rriQ, and a characteristic length, l/rriQ, will be used to obtain di- 

mensionless quantities in what follows. By considering a certain path partition {zq, zi, , zat-i, ^at, z), 

we can write: 

p e'/; Qir) ^ ^ ^11, dr QM^ ^^ ^^^^ dr QM ^ /^^ Qir) ^^^^ 

As shown in the appendix, we find for small Aj, where Aj = Zj^i — zj, that 
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P = e^-3P.^. + 1 ^ J. sin(p,A,) . (13) 

where pj = +{E'^ — fjY^'^ and fj is defined as 

f^ = ^J^l'''' drfir), (14) 

e.g. as the average value of the integraL 

Thus, by substituting the path order products of the right-hand side in (|T^ by (|TB|), 
multiplying and reordering the terms, we obtain 



n(z, zq) = e + 

N k-1 

+ -a2j2MPk^k)fke ^='=+1 e +0 

^ k=0 



1' 
E 



(15) 



where we have used e^A = Ae~^ when {A, B} = 0. If the path partition is now considered 
to be thiner and thiner, or in other words, if the limit ^ is taken for all j, then 
we can replace ~^ S dr. The definition of fj as an integral average value leads to 



the following substitutions fj f{r) and pj — > p{r), where p{r) = + [e'^ — [/(r)]' 
provided that /(r) is continuous. In this way, we obtain 



1/2 



Vl{z, Zo) 



^n, I drp{r)f{r)e-^"''rr<^PiO 



+ 



1' 
E 



o J zn 



(16) 



We must here stress that VL{z,zq) has been expanded in powers of 1/-E, but not in a 
strict sens because the coefficients are also depending on E. As a consequence of this fact, 
whether the series for VL{z, zq) which is truncated in order to give (|16]) is an asymptotic 
one or not must be carefully elucidated. This point is treated in detail in |§] and we will 
return to this question below. If we define the quantity a = 2aE and expand all the terms 
which depend on E in ( |16|) in powers of 1/E, by assuming a constant, we obtain 



il(Xa, —Xa) 



1 aX 

1 + r^T 



dx F{Xx) e 



-icr^aXil — x) 



+ 



E 



(17) 



where z = Xa and zq = —Xa, with A ^ 1. The function F{r) has been introduced in (|TT|) 
and the integration variable is changed to x = j^. The quantity is defined as "classical 

path" , drpij) = Xa + O 



(;!) . Equation ([TTD can be written in a matrix 



way 
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n{Xa, -Xa) = ( ] . 

Neglecting terms of the order o{l/ E) we get 



LOll — " — 1^22 5 

= — ^^^w = ^2i; (19) 

where 

1(a) = y^^dxe*^"^F(Aa;). (20) 

By considering the requirements imposed to the profile wall function, /(r), the follow- 
ing approximations are performed: m(r) =0 for t < —Xa and m(r) = mo for r > Aa; i.e. 
F{Xx) = 0, 1 for X < —1 and x > 1, respectively. Thus, taking into account the chirality 
eigenvalues of the eigenspinors tpi, ijj2, V^s and ip^, we have that tpi and ^/'2 correspond 
to right-moving right-handed particles and left-moving left-handed particles respectively, 
and ips and 7/^4 to right-moving left-handed and left-moving right-handed particles, in the 
region where x < —1. Concerning the region where x > 1, Q{Xax, Xax) can be im- 
mediately diagonalized in order to identify the right-moving and the left-moving flux of 
particles, since the mass is constant. Thus, by requiring that the left-moving flux is zero 
in the latter region, i.e. by imposing that we only have a transmitted flux in this region, 
a left-handed reflected flux is obtained from the right-handed one and vice versaQ. The 
equivalence between the equations for ipj and for ipn in (^) , as a consequence of the result 
Q{z) = Q{z) for the real mass case, implies that there is no asymmetry in the reflection of 
incident right-handed and left-handed fermions. In consequence, needless to distinguish 
ipi and Only when an imaginary mass term, which introduces CP- violating effects, is 
incorporated, differences between the behaviour of right-handed and left-handed fermions 
may exist in the reflection process[l,3,5]. The fact that CP invariance boils down to an 
equality between both helicities is a consequence of CPT invariance [0. 

In this way, following the method above described, which was originally introduced by 
Nelson et al for a numerical computation of the reflection coefficient [6,7], we use that 

tlj{z)=n{z,ZQ)ij{zo), (21) 

where 



H^) = D[ I and iji-z) = . (22) 



^Indeed, by conservation of angular momentum the helicity is flipped in the reflection process. 
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D diagonalizes the matrix Q in the region where the mass is constant. From eq. 
fB^) we obtain for the reflection coefficient 



uli - {E - p)Ui2 



with p = +\/E'^ — 1. In our case, for the high energy limit, we have 



R{E) = ^ (e'^'^ - z\al{a)) + o 



(24) 



Now, we must compute the integral which gives the function X(a) in eq. (pOl). We 
introduce a certain convergence factor, e~^^, where the parameter e is taken small enough 
to be sure that E{Xx) decreases faster than this factor when x — oo, in such a way that 
the integration over the real axis, from — oo to oo, is well defined. Thus, by cutting the 
integration domain into three intervals, (— oo, —1), (—1, 1) and (1, oo), we can write 

J(a) = lim Qy°° - - y ^|cixe(*^"~")^'F(Aa;)^ . (25) 

where the integrals in the regions (— oo, —1) and (1, +oo) can be trivially solved. Indeed, 
as indicated above, we approximate E{\x) = 0, 1 for x < — 1 and x > 1, respectively. 
Therefore we have 



lim 



lim 

£^0+ J + 1 



-OO 
+ 00 



dx E{Xx) e(*^"-^)" = 
dx F{Xx) e(*^"-")^ 



Xa 



(26) 



Concerning the integral over all the real axis, the Cauchy theorem can be applied in 
order to perform this integration. If the Laurent expansion for F{z) in the pole of order 



Zj, is Hii.^ ~ -^i)" the following result can be obtained 



lim 



+ 00 



dx F{Xx) e 



{i\a—e)x 



2Tii 
~X 



N 



[la 



,n-l 



(n-l] 



(27) 



j=l n=l 

where zi = xi + iyi, Z2 = X2+iy2, zn = XN + iyN, are all the poles of F{z) with positive 
imaginary part we have picked when the integration contour is adequately closed. Since 
the profile function is analytic on the real axis, it is obvious that yj ^ and therefore we 
obtain an exponential dependence on a. 

It can be shown from eqs. (p4|), (^), (|^) and (^71) that 



■jrf] 

R{E) = — V e-^y^ e'"''^ 

^ 3=1 



N 



[la 



in— 1 



n=l 



(n-1)! 



+ 



e) 



(28) 



Obviously, R{E) defined in eq. ( PHj ) and given by eq. (^) is the reflection coefficient 
generated by a profile function which is zero in the region (— oo, —Act) and constant in the 
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region (Act, oo). Nevertheless, /(r) is an analytic function in the real axis verifying that 
/(r) — 0(t) decreases exponentially when r — > ±cx), in such a way that the wall profile is 
well defined over all the real axis. In fact, we are interested in the situation where the wall 
profile extends from —oo to oo, although the former requirement is technically necessary 
in order to identify the incident, reflected and transmitted waves. On the other hand, we 
cannot take directly the limit when A ^ oo||^. This problem is solved by considering A 
as large enough to neglect the effects due to the evolution of the wave outside the region 
(— A(T, Aa). The question is studied in detail way in ref. [H and the result is, as intuitively 
expected, that in the limit where we take for the wall profile an analytic function from 
— oo to +00 one simply has to consider in ( ^5|) the first integral, from —00 to +00, and 
that consequently p8| ) gives the reflection coefflcient for the wall from —00 to +00. The 
effect of the mentioned technical assumption (the substitution of the true profile function 
by a step function outside the interval), is that T(a) is defined as an integral over de 
domain (—1, 1). It can be bounded as seen in (^5]) when this function, X(a), is expressed 
through the integral over all the real axis and the two other integrals given by eq. (PBj). 
It is worth to stress that indeed R{E) in eq. (p8| ) does not depend on A, as expected for 
the reflection coefficient for the wall from —00 to +00. 

In (^) the leading term for a series of powers in has been written. Nevertheless, 
we can be sure that it gives the leading term only if a is a constant. If the parameter 
which is considered constant is a, then we must return to the problem of the asymptotic 
properties of the series we use. As above mentioned, this question is discussed in 0, 
where we conclude that the right-hand side of (^) is the leading term of the reflection 
coefficient for high energy when o" -C 1, even if a = 2(7 E S> 1. In this way, by considering 
the dependence on E for the reflection coefficient in the high energy limit, taking a as 
constant, we have 

R{E) = 2'Ka J2 e-^'^"''^ e^^^^^^^ ^-"^T^^^ • (29) 

j=l n=l -'-)• 

Thus, the behaviour of the high energy fermions, i.e. the high energy component of a 
wave packet, hitting the wall is completely determined by the complex singularities of 
the wall proflle function. Furthermore, if we are in such a range of energy that aE 3> 1, 
the leading contribution is exponentially decreasing with the energy, the coefficient of the 
exponential being given by the pole of f{z) with the smallest imaginary part, i.e. the 
closest to the real axis. Calling = Xk + iyk this prevailing pole we obtain 



{2iEa) 



R{E) = 27rae~2^'^^^ bt,^ ^ ^ e^^^^^^^ . (30) 

The dependence in a in (|30|) is simply an artifact of our change of variables (|l3)- Indeed 

+00 , . 

it can be inmediately seen that if we consider the Laurent expansion, b^{z — Zj)'^, 

n=—u 

for f{z) instead F{z), in the pole of the order uj, z = Zj, we have 
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where, by using ([TI|) , the poles and the coefficients of the Laurent expansion for F{z) and 
f{z) have been related. Thus, the final result can be expressed in terms of the singularities 
of f{z) by replacing (|3ll) in eq. and (pOD, 



.=1 



where the lower expression is valid in such a rang of energy that e"^'^*^^*"^''^ can be 
neglected, and being the imaginary parts of the two closer poles to the real axis (this 
requirement is analogous to the former aE ^1). The dependence on a has disappeared 
in (|32|) , as expected. 

The results we present can be checked by using the particular analytic solution ob- 
tained in ref. |^, Q for the following wall profile 

/(r) = i±i5!^ . (33) 



For this particular proffie function, it follows from ( P9|) by computing the sum for the 
infinite series of poles it presents 

R(E) = , (34) 

^ ^ 2sinh(7rEa) ' ^ ^ 

which agrees with the result of reference |0 and in the high energy limit for a ^ 1. 

With this positive check, we conclude that (|29|) in general and (|30D for Ea large give 
the high energy behaviour of the reflection coefficient for analytic profile functions in the 
real axis. The singularities of the profile function determine the behaviour of the high 
energy fermions in the way shown by eq. (pOD , and if the energy is larger enough than 
the inverse of the characteristic wall thickness, it is the closest pole to the real axis that 
gives the high energy behaviour through eq. (|30D . 
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A Appendix 



We take Q{t) as constant in the intervals {zj, -Zj+i), defining for eacli one fj as tlie following 
integral average value 



1 /"^j+i 



(35) 



where = 2:^+1— 2;^-. By assuming that Aj is small (in fact, the limit — is considered 
for the result obtained in this appendix), it can be written in good approximation by using 

(I) 



i n + l dr Q(r) r 



Taking into account that for two operators verifying {A, B} = and = B'^ = 1, it can 
be easily proven the following result 



nA+l3B 



cosh 



1/2 



aA + pB 
H 77? smh 

(^2+ ^2)1/2 



{a' + 



1/2 



(37) 



and from (|36| 



P ^^nr "^^^^ = e-3P.^. + 1 sin(p,A,) + O 

E 



e) 



.(13) 



where pj = +{E'^ — fjY^'^. In order to obtain ([T3|) ^ and have been expanded in 
powers of ^. 
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